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We consider a quantum mechanical particle living on a graph and discuss the 
behaviour of its wavefunction at graph vertices. In addition to the standard (or 
S type) boundary conditions with continuous wavefunctions, we investigate 
two types of a singular coupling which are analogous to the 5' interaction and 
its symmetrized version for particle on a line. We show that these couplings 
can be used to model graph superlattices in which point junctions are replaced 
by complicated geometric scatterers. We also discuss the band spectra for 
rectangular lattices with the mentioned couplings. We show that they roughly 
correspond to their Kronig-Penney analogues: the 5' lattices have bands 
whose widths are asymptotically bounded and do not approach zero, while 
the 5 lattice gap widths are bounded. However, if the lattice-spacing ratio 
is an irrational number badly approximable by rationals, and the 5 coupling 
constant is small enough, the 5 lattice has no gaps above the threshold of the 
spectrum. On the other hand, infinitely many gaps emerge above a critical 
value of the coupling constant; for almost all ratios this value is zero. 



1 Introduction 

The problem of describing a quantum particle living on a graph is not new in quan- 
tum mechanics; it appeared for the first time in early fifties in connection with the 
free-electron model of organic molecules p6 |. Writing down a Hamiltonian of such a 



particle requires to check that the coupling between the wavefunctions at branching 
point of the graph defines a self-adjoint operator, or in physical terms, that the 
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probability flow is preserved there. This is conventionally achieved by demanding 
that the wavefunctions are continuous at the junctions and satisfy the conditions 



E 
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f'AXm) = C m f(x m ) , (1.1) 



where m is the vertex number, the sum runs over all links entering this vertex, 
f(x m ) is the common value of the functions fj there, and real parameter. 

In recent years the interest to quantum mechanics on graphs has been revived 
- see, e.g., Jj], || ||, || || [R| [19|, ^(J and references therein — in particular, as a 



reaction to the rapid progress of fabrication techniques which allow us nowadays to 
produce plenty of graph-like structures of a pure semiconductor material, for which 
graph Hamiltonians represent a natural model. This posed anew the question about 
physical plausibility of the boundary conditions (|1 . 1|) . 

This problem has two basic aspects. The first of them concerns the derivation 
of graph models from a more realistic description, in which the configuration space 
consists of a system of coupled tubes. This is still not the true system in which 
the tubes are complicated many-body objects, but the crystallic structure of the 
semiconductor material makes it a reasonable approximation. 

Replacing a tube system by a graph of the same geometry, say, by the tube 
axes, means a substantial simplification. The reason is that one cannot assume that 
the wavefunction in a tube is independent of the radial a azimuthal coordinates as 
the authors of |HJ naively suggest. The tube Laplacian is specified by appropriate 
boundary conditions, so there is a fully concrete system of transverse eigenfunctions; 
most frequently the tube boundary is assumed to be Dirichlet in which case none of 
these eigenfunctions is constant. The graph approximation is generally expected to 
work in situation where the tubes under consideration are thin enough so that the 
transverse-mode eigenstates are well separated in energy and their geometrically 
induced coupling coming from bending and branching is weak. 

The one-mode approximation for a single bent tube can be justified [T^, [L3| but 



the problem is in no case a simple one. The case of a branched tube is even harder 
because a typical branching region is self-similar with respect to changes of the tube 
widths, so there is no natural parameter to switch off the intermode coupling. In 
general, we know neither the condition under which the graph approximation works 
for a branched tube, nor the boundary conditions which should model a particular 
branching geometry. However, we are not going to discuss these problems in the 
present paper, apart from some comments given in the concluding section. 

The second aspect concerns intrinsic properties of the graph Hamiltonians: one 
may ask what is the family of admissible operators and which place is occupied in 
this class by those referring to the conditions ( |1.1|) . This problem was solved in |16| 



where we showed how a general graph Hamiltonian can be constructed using the von 
Neumann theory of self-adjoint extensions. However, the operator family obtained 
in this way is large: even if we exclude nonlocal interactions, i.e., we forbid particle 
hopping between different branching points, each vertex of the graph is associated 
with n 2 real parameters, where n is the number of links entering the junction. 
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Motivated by this we discussed in |L6| in detail several subfamilies of such opera- 



tors. The simplest situation occurs if the domain of a graph Hamiltonian is required 
to consists of functions which are continuous at each vertex. Then we arrive back at 
the situation we started with, namely at the conditions ( |1 . 1|) in which every junction 



is characterized by a single parameter. Contrary to the claim made in [19], however, 
this parameter is a real number and not a function of energy. 

The above discussion shows that until derived within a squeezing-tube approxi- 
mation, the wavefunction continuity is just an assumption which we may or may not 
adopt. If we decide to drop it, the next more general class of graph Hamiltonians 
consists of those which are locally permutation invariant at each vertex; we have 



shown in [IB] that any junction is then described by a pair of real parameters in 
such a way that 

fj( x m) = A m fj(x m ) + B m ^ f'k( x m) 3 3 = ^-i---i n mi 

where the indices have the same meaning as above; this condition reduces to ( |1 . 1| ) 
for A m = B m =: c" 1 . At the same time, a junction can be described also by two 
one-parameter families of boundary conditions, which represent singular limits of 
(|1.2|) when A m , B m — > ±oo while 

C m := B m — A m or D m := n m [A m + (n m — l)B m ] (1-3) 



is kept preserved; they are described by the boundary conditions ( |2.5| ) and fl2.6|) 
below, respectively. 

Although these couplings and their scattering properties were discussed in detail 
TE \, their meaning remained somewhat obscurred. Our aim in this paper is to 



in 



show that they represent in a sense a counterpart to the coupling given by (|1 . 1| ) 
and that they generalize naturally the concept of 5' interaction which has attracted 
attention recently in connection with spectral properties of Wannier-Stark systems 
H |7|, ^3|. This will be done in the next section; we are also going to show there 
that the relation between 8' and geometric scatterers discovered in H] extends to 
vertices of any number of links. 

After that we shall discuss the spectra of periodic rectangular lattices in which 
each junction is described by one of the couplings mentioned above. For the sake 
of simplicity, we restrict ourselves to the planar case, however, the results have a 
straightforward extension to higher dimensions. We shall demonstrate, in particular, 
that the 5 lattice spectrum depends substantially on the ratio 9 of the rectangle 
sides. It can even happen that there are no gaps above the threshold of the spectrum; 
this is the case if 9 is badly approximable by rationals and the coupling is "weak" 
enough. On the other hand, infinitely many gaps exist for almost all 9 but the 
gap pattern, as well as the band pattern for a S' s lattice, is again irregular for an 



irrational 9 . A summary of the results has been given in |TJ . 
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2 5' and 5' s interactions 



Throughout this section we will consider a single branching, hence we may drop the 
index denoting the junction. For simplicity we will deal with the graph T n consisting 



of n halflines whose endpoints are connected at a single point; as explained in |L6 
the boundary conditions describing the coupling are local and may be used for any 
n-link vertex. 

The state Hilbert space for a spinless particle having T n as the configuration 
space is H := 0" = i L 2 (0, oo) and the free Hamiltonian acts at / = {fj} as Hf = 
{—fj} ', its domain consists of all functions whose components are H 2,2 on the 
halflines and satisfy appropriate boundary conditions at the connection point which 
relate the boundary values fj(0+) and /j(0+) , j — 1, . . . , n . For the sake of brevity 
we shall drop the arguments as long as a single connection point is concerned. 

2.1 A warm— up: singular interactions on line 

Consider first the case of T2 which is naturally isomorphic a line with a single point 
interaction. In general, there is a four-parameter family of such interactions — cf. 
JJ, [K], [U], |29|. The best known among them are the 5 interactions for which 



the wavefunction is continuous at the point supporting the interaction, while its 
derivative has a jump proportional to the function value 0, Sec. 1. 3]. Another im- 
portant class are the 8' interactions where the roles of the functions and derivatives 
are reversed: the wavefunctions now satisfy the conditions 

/; = /!:=/', f+-f- = pf, (3 em. (2.1) 

The name is somewhat misleading because in distinction to 8 , the 8' interaction 
cannot be approximated by Schrodinger operators with scaled potentials [||, |30|. In- 
stead, one can use approximations by families of rank-one operators |3(J or velocity- 



dependent potentials [IU|. Another way to understand the 8' interactions has been 
suggested in 0: replacing the line by a "spiked-onion" graph, i.e., cutting it into 
two halflines and joining the loose ends by N links of length I , one can reproduce 
the high-energy scattering behaviour of the 5' interaction up to a phase factor in 
the limit when N — ■> oo , i — > and (3 := Nl is kept preserved. 

In the two-parameter family of space-reversal invariant interactions there is 
another subset with a similar property which may be regarded as a symmetrized 
version of the 5' interaction; it is characterized by the boundary conditions 

f + + f- = 0, f + + f_ = Df + , DER. (2.2) 

It is easy to check that it has a bound state of energy — 4D~ 2 provided D < and 
the reflection and transmission amplitudes for a plane wave of momentum k are 

2 — ikD 2 — ikD 

which up to signs coincides with both the 5' result (for (3 = D) and the high-energy 
behaviour of the limiting geometric scatterer. 
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2.2 Vertices with n links 



Let us look now how the above results extend to a quantum motion on T n with 
n > 3 . The coupling with continuous wavefunctions is a natural analogue of the 5 
interaction, so we shall use for it the same name. The boundary conditions ( |1 . 1|) 
can be rewritten as 



n 

/! = ••• = /»=:/, E/; = c /- (2-4) 

3=1 



The first of the two classes mentioned in the introduction is for a given C G M 
characterized by the relations 

n 

E/j = and f3-fk + C(f' j -&) = 0, j,k = l,...,n; (2.5) 

i=i 

it is clear that just n among these conditions are independent. For n = 2 this 
coincides with the requirements Q2.1| ) if we put (3 = —2C ; recall that both links are 
now positive real halflines so one of the derivatives has to change sign. In view of 
this analogy we shall refer to the coupling ($.dj as to 5' . It is sufficient to consider 
C / 0, because otherwise the conditions fl2.5|) reduce to ( |2.4j ) for c = . On the 
other hand, the second exceptional class is given by 

n 

/[ = ■■■ = fn =:/', Efi = D f ( 2 - 6 ) 

i=i 

with D G iR . In the two-link case this correspond to the condition ( |2.2|) , hence we 
shall speak about the 5' s coupling. 

To justify the claim represented by the relation (|1.3j ), let us solve (|L2 ) with 
respect to the derivatives. This yields 

f] = (A -B)V n -i fj ~ (A-B)V n ^ g. (2 ' 7) 
where T> m := A + mB . In the first limit of ( |1.3| ) we arrive then at 

cf; = —fj + -E/*> 

J n n rr. 

which is equivalent to fl2.5Q, while the other limit gives 

1 n 

fj = p E /*» j = 1, . . . , n , 
fc=i 

i.e., the condition Q2.6Q . 

Each of the three operators on r„ described above can exhibit bound states with 
eigenf unctions localized around the junction — this happens iff the corresponding 
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parameter, c, — C or D , respectively, is negative (we keep the notation of [|T(| ; it 
would be more natural in the present context to change the sign of C). Since any 
solution to the free Schrodinger equation is of the form {aj e~ KXj } , it is sufficient to 
substitute to the boundary conditions ( ^4|) -( ^6|) to check that the corresponding 
eigenvalues are 

6(c) = - (£f , e(C) = -C- 2 and e(D) = - 

in the three cases; the 6 and 5' s eigenvalues are simple, while the 5' bound state 
has multiplicity n — 1. Also the scattering properties of these junctions are easily 
found — see [16| . 



2.3 Geometric— scatterer junctions 

Now we want to show that the "spiked-onion" argument mentioned in Section 2.1 
can be naturally extended to 5' and 5' s couplings on T n , n > 3 . Let us replace the 
latter by the graph T n (N,£) sketched on Figure 1: every pair of hamine endpoints 
is connected by N links of length i ; the corresponding variables will run through 
the interval [— £/2, £/2] . We shall assume that the coupling at each graph node is 
given by the condition ( |2.4| ) with the same parameter c . 

Using the permutation symmetry of the Hamiltonian, we can write a general 
scattering-solution Ansatz in the following form: 

(i) e~ lkx 4- re ikx a t a chosen external link (hamine), 

(ii) te lkx at the other n — 1 halffines, 

(iii) a e lkx + (3 e~ lkx at the 2N connecting links coupled to the "incident" halffine, 

(iv) 5 cos kx at the remaining iV ^) ~ 2 connecting links. 
The condition ( ^.4| ) now yields the following system of equations 

1 + r = afj + f3rj 



t = arj + (3rj = 5 cos 
r — 1 + N (n—1) (afj — (3rj) = 7(1 + r) 



hi 

t- N(arj- /3fj) -iN(n- 2) 5 sin— = 7*. 



(2- 



where we have denoted 



c 
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The first three equations coming from the continuity requirement are solved by 

1 fj{l+r) — rjt i fjt — 77(1+ r) ^ t 

2 sin k£ 2 sin k£ cos y 

Next we substitute these values into the remaining two of the equations (|2.8|) and 
use the identity 

cot 2x — (n— 2) tana; = (n— 1) cot 2x — (n— 2) esc 2x ; 

introducing 

P = p(k) := 1 - 7 + ziV(n-l)cot£;£, Q = Q(k) :-- 



sin ki 

we can rewrite the resulting system as 

rP - (n-l)Qt = P , -rQ + t[P-(n-2)Q] = Q. 
Since Re P = 1 , the sought reflection and transmission amplitudes are given by 

|P| 2 - (n-2)PQ + (n-l)Q 2 



r(k) 
t(k) 



P 2 - (n-2)PQ + (n-l)Q 2 ' 

(2.9) 

2Q 



P 2 - (n-2)PQ + (n-l)Q 2 ' 
It is straightforward to check the identity 

\P 2 - (n-2)PQ + (n-l)Q 2 \ 2 - \\P\ 2 - (n-2)PQ + (n-l)Q 2 \ 2 = 4(n-l)|Q| 2 , 
so the S-matrix is unitary, 

|r(A;)| 2 + (72-l)|t(A;)| 2 = 1. (2.10) 

Consider now the limit of increasingly complicated scatterers; in other words, we put 
£ :— and let iV — > 00 . The expressions ( |2.9| ) then can be for large N rewritten 



as 



n -2 - S£ - (f)ik T 



(2.11) 



from where the limits immediately follow. At large energies, in particular, the terms 
containing c can be neglected and the amplitudes ( |2.11| ) behave as 



n-2 - [ n AikT _ 2 
r(k) « j^l— , » j-— ; (2.12) 



-n 



+ (fjikr ' ' -n+ (£)ikr 



if c = , of course, the right sides give the expressions of the limiting amplitudes 
for any k . 

Let us now compare this result to the S-matrices corresponding to our two 



singular couplings which have been computed in [16]. In the 8' case the reflection 
and transmission amplitudes are 

2 — n + inkC .,. 2 . 
r ( k ) = ,^ . t( k ) = 2 - 13 

choosing r := — , we get at high energies the same t(k) , while the reflection 
amplitude differs by the phase factor —2 arg(n — 2 + inkC) which goes to tc as 
k — > oo . As for the 5' s coupling, we have 

/,x n — 2 — ikD —2 

= T7^' '(* = TH' ( 2 - 14 ) 

n — z/cD n — iku 



so one has to put r := n ^_^ to recover the S-matrix elements (|2.12| ) up to a sign, 
which is switched for all of them. 

Hence both the singular couplings reproduce the high-energy behaviour of the 
limiting geometric scatterer up to a phase factor. They are the only ones with this 
property in the class ( |1.2|) . To see this, consider the reflection amplitude correspond- 
ing to ( |1.2j ). Leaving aside the case A = B corresponding to fl2.4| ) when we have 
lim^oo r(k) = , its high-energy behaviour is 

n - 2 + ik (l - 4) Pft-i 

r(k) = , + 0(k~ 2 ) , 

-n + 2(l-i)+ik(l-i)v n „ l 



which certainly differs from ( |2.12| ) by more than a phase factor. On the other hand, 



all these couplings (with exception of the case A — B ) represent the effective Neu- 
mann decoupling at high energies, lim^oo r(k) = 1 , while the geometric scatterer 
mimicks rather the Dirichlet decoupling, limfc^ 00 r(fc) = —1 . 



3 Lattices with a singular coupling 

Consider now the case which has attracted some attention recently [0, [l9j as a 
model of quantum-wire superlattices. We shall assume that the graph in question is 
a rectangular lattice with the spacings £\, £2 in the x and y direction, respectively 
(cf. Figure 2). In addition, we suppose that each graph vertex is endowed with 
the same coupling of one of the above described types. We restrict ourselves to 
the planar situation just for the sake of simplicity; the band conditions obtained 
below and the method of their solution have a straightforward extension to higher 
dimensions. 



8 



3.1 The Bloch analysis 

To find the band spectrum of such a lattice, we start from a natural Ansatz for the 
Bloch solutions: we choose 



f m (x) = e m9 ^ (a n e lkx + b n e~ tkx ) ... x G (n£ h (n + l)h) 

9n(y) = e iM (c m e^ + d m e- lk y) ... ye {m£ 2 , (m + l)£ 2 ) 
The coefficients on neighbouring links are related by 

a n+1 te ikx + b n+1 Ze- ik * = a{a n e ik * + b n e~ ikx ) , 

c m+lV e ik y + d m+1 fje- ik y = r(b m e ik y + d m e~ ik y) , 

where we have denoted for the sake of brevity 

a ■= e i9l£l , t := e ie2h , f := e ik£l , 77 := e ikh ; 

this conditions allow us to compute easily the needed boundary values. 
Let us begin with the 5 coupling, where the relations fl2.4p now read 



(3.1) 



(3.2) 



fmirdt+0) = f m {n£i-0) = g n {m£ 2 + 0) = g n {n£ 2 -0) =: F mn , 



4(^+0) - + g' n (m£ 2 +0) - g' n (n£ 2 -0) = cF n 



(3.3) 



Substituting the boundary values, we get a homogeneous system a four independent 
equations for the coefficients, which has a solution provided 

-2iA;f(?7-77)[l -(!(£+£) +<x 2 ] + ccrf(£- 0(77-77) - 2ika^-0[l-f(7]+f]) + f 2 } = 0. 
Returning to the original quantities, we can cast it into the form 

cos^i^i — cos Hi ^ cos 9 2 £ 2 — cos k£ 2 c Q- (34) 
sin k£i sin k£ 2 2k 

this is the result obtained in |TD[ , and in H for the particular case c = . 



If we assume instead that the lattice links are coupled by the 5' s interaction, the 
requirement (|3~^ ) is replaced by 



run j 



(3.5) 



fL(nii+0) = -f'Jrdx-ti) = g' n (m£ 2 + 0) = -g r n (n£ 2 -0) =: G r> 

+ fmin^-O) + g n (m£ 2 +Q) + g n {n£ 2 -0) = DG 
solving it in the same way we arrive at the condition 

cos^i^i + cos/c^i cos 8 2 £ 2 + cos k£ 2 Dk 

sin kit + sin kl 2 2~ = ' ( ' ' 



Finally, in the 5' case the coupling conditions read 

f^(n£ 1 +0)-fl n (rd 1 -0)+g' n (m£ 2 +0)-g' n (m£ 2 -0) = 0, 

fU^+0)-g n (m£ 2 +0) + C(fl n {ne 1 +0)-g! n (m£ 2 +0)) = 0, 

(3.7) 

g n (m£ 2 + 0)-f m (n£i-0) + C(g' n (m£ 2 + 0) + f^n£ 1 -0)) = 0, 

/mK-o) - ffn (mi 2 -o) + c(-/:k-o) + 9 ;k-o)) = 0. 

The solvability condition is now slightly more complicated; it can be written as 



' cos^--Re((l-ifcC)e^) 
^ Im((l+^C , )-V fc£ ^ " ' ' 



At high energies, however, this relation simplifies: up to 0(k 1 ) terms it acquires 
the form 

cos^i^i — cos Hi cos 9 2 £ 2 — cos k£ 2 . 

— + ±m (M) 

For a single junction considered in Section 2.3, a 5^ coupling with the parameter 
D has, up to a phase factor, the same asymptotic behaviour as 8' for C = —D / n . 



Using this substitution for n = 4 we arrive at the relation (|3.6| ) with the parameters 
9j£j replaced by 8j£j + tt . Hence the band spectra of the 5' and 5' s lattices behave 
alike at high energies; we shall discuss below only the second case which is simpler. 

3.2 S lattice spectra 



The condition ( |3.4j ) yields no restriction on k if c = as the authors of || pointed 
out; the spectrum covers the interval [0, oo) and in the isotropic case, £\ — £ 2 , one 
can write the energy e(k 2 ) := k 2 in terms of the Bloch parameters (quasimomentum 
components) 9i, 9 2 explicitly. 

This is no longer true if the coupling constant is nonzero. Nevertheless, one can 
say a lot about the spectrum determined by ( |3.4|) . To solve this condition, let us 
rewrite it as 

^(M-x;^^, (3.io) 

3=1 J 

where Vj := cos9j£j . Since these parameters run through the interval [— 1, 1] , we 
find easily 

k£j Vi cos k£i k£j , „ „ 

- cot < - 1 — - 1 < tan — £ ... sin kL > , 

2 ~ sinA;^ 2 3 

Ma Un cos k£-j k£j 
tan — 1 < 3 . , - 3 < - cot — 3 - ... sin kL > . 
2 ~ sinfc^- ~ 2 3 
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From here we obtain the extremum values of F(k; •, •) for a fixed k : 



F + (k) := max F(k;vi,v 2 ) 

iye[-l,l] 



F^(k) : = min F(k;vx,V2) 



tan 



(Mi 


71 


k£j 




{- 


~~ 2 


TX 


) 



(3.11) 



(kL 


71 


k£j 




{- 


~ 2 


71 


) 



where the symbol [•] denotes the integer part; at the singular points the value is 
understood as the limit from the right or left, respectively. By definition we have 
±F±(k) > , hence the gap condition can be expressed as 



± Tk > ±F ± (*) 



(3.12) 



for ±c > , respectively. This applies to positive enegies, on the negative halfline 
one compares instead j- , where k := — ik , with the extremum values 



F + (k) := - tanh 
j=i 



txlj 



F_(k) := - E coth 



(3.13) 



since both of them are negative, there is obviously no negative spectrum for c > . 
Notice that the band condition for a one-dimensional array of 5 interactions can be 
rewritten in the same form with a single term on the right sides of the relations ( ft.lip 
and (|3.13|) ; it is easy to reproduce from here the spectrum of the Kronig-Penney 
model §, Sec.III.2]. 

Let us collect first some simple properties of the spectrum which follow directly 
from the condition fl3.12j ) and its negative-energy counterpart. We introduce 



so 



•1*2 ) 



-1/2 



e :-- 



and L := max(£ l5 £ 2 ) 



(3.14) 



£9 1 ' 2 



min^A 1 ) 



L- 1 



and denote by a(£, 6, c) the spectrum of the corresponding 5 lattice Hamiltonian. 

Proposition 3.1 (a) The spectrum has a band structure, a(£,9,c) = |j£Li 
for some N > 1 , where a r < (3 r < a r+ i . 

(b) a(£,9,0) = [0,oo) . 

(c) If c > , each f3 r is of the form (^f) f or ^ = £i or £ 2 , and some n G /Z . 
Similarly, a r = (j^) with m G for c < and r > 2 . 

(d) ±«i > holds iff ±c> . 
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(e) For c < -4r x (0 1/2 + O^ 2 ) , fa < and a 2 = (if . 

(f) a{£,6,c')nlR + C a{£, 9, c) n R+ holds for \c'\ > \c\ . 

(g) Each gap is contained in the intersection of a pair of gaps of the Kronig-Penney 
model with the coupling constant c and spacings t\ and £2 ■ 

(h) All gaps above the threshold are finite. If there is an infinite number of them, 
their widths are asymptotically bounded, 

a r+1 - (3 r < 2\c\{£ l +£ 2 y l + 0( r - 1 ) . (3.15) 

The property (g) allows us to localize spectral gaps by those of the Kronig- 
Penney model. By negation, it illustrates that transport properties of the lattice 
are better than a combination of its x and y projections. If a given energy is 
contained in a band in one of the directions, then by ( |3.4| ) it is trivially also in a 
band of the lattice Hamiltonian. The converse is not true, of course: the condition 
may be satisfied even if none of the factors can be annulated separately. Of course, 
different solutions correspond to different quasimomenta and different directions in 
which the particle is able to "dribble" through the lattice. 

However, the most interesting property of the spectrum is its irregular depen- 
dence on 9 coming from the existence of competing periods in ( [3.11|) . To formulate 
the results, we have to recall some notions from the number theory [[H], |2"T| . Irra- 
tional numbers can be classified with respect to how "fast" they can be approximated 
by rationals. In particular, such a number is called badly approximable if there is a 
5 > such that 

5 

> 3.16 

q q 2 

This is a non-empty subset in the family of all Diophantine numbers; for instance, 
it contains all algebraic numbers of degree 2 , i.e., irrational solutions of a quadratic 
equation with rational coefficients. On the other hand, the Lebesgue measure of this 
set is zero. 

One can also write 9 as a continued fraction [ao,oi,...] with integer coeffi- 
cients; such a representation is unique and provides a natural way to gauge the 
approximation properties. The faster the a n 's grow, the better is 9 approximated 
by the truncated fractions; the worst irrational from this point of view is the golden 
mean | (1 + v^5) = [1, 1,1,...]. A counterpart to badly approximable numbers is 
the class of irrationals with an unbounded sequence of coefficients, 

limsup a n = 00 , 

n— >oo 

which has the following property: there are sequences {m r } ( ^L l , {n r }'^L l of pairwise 
relatively prime integers such that 



9- P - 



lim n r 

r—>oo 



9 - ^ 

n r 



0. (3.17) 
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These numbers, which may be called Last admissible p2 |, have full Lebesgue mea- 
sure. It is sometimes convenient to express these approximations using number- 
theory symbols: the fractional part {x} := x — [x] , and ||x|| := min({x}, 1 — {x}) . 
Since the number in question may be indexed, to distiguish the fractional part from 
a sequence, we shall always specify the range of indices in the latter case. 

Theorem 3.2 (a) If 6 is rational, o~(£,8,c) has infinitely many gaps for any 
nonzero coupling constant c . 

(b) For badly a approximable 9 there is cq > such that for |c| < cq the spectrum 
has no gaps above the threshold, Pi — oo . 

(c) a(£,9,c) has infinitely many gaps for any 9 provided \c\L > 5 _1 / 2 7r 2 . 

(d) If 9 is Last admissible, there are infinitely many gaps for any c^O. 



Proof: We shall suppose that c > ; the argument for a negative coupling constant 
is analogous. If 9 — | , the function F + has infinitely many zeros in M + without 
accumulation points, so (a) follows. On the other hand, F + (k) > for any k > 
if 9 is irrational, hence we have to investigate local minima of this function. They 
occur at k n := and k m : = ^ with n,m <G ZZ + , and the corresponding values 
are 

F + (k n ) = tan(|{n0}) , F + (k m ) = tan (^{rnt 1 } 



If 9 is badly approximable, the condition ( |3.16|) yields 



F+(k n ) > -{n9} > -\\n9\\ > — . 

+ 1 n) 2 \ { - 2 || II ^ 

Since 9~ x is also badly approximable, we get F + (k m ) > ; if different constants 
correspond to 9, 9" 1 , we call 5 the smaller of the two. This implies that for c small 
enough 



2k n > -^' m ' 2k n 
holds for all n,m G /Z+ , i.e., the assertion (b). 



By the Hurwitz extension of Dirichlet's theorem [27, Sec. II. 1] one can find to any 
irrational 9 sequences {n r }^ 1 , {m r }'^L 1 in such a way that \n r 9 — m r \ < h~ l l 2 n~ l . 
Moreover, these approximations can be constructed explicitly in terms of truncated 
continued fractions |2T|, Chap. 10]. Choosing the truncations of even lengths (without 
the integer part), we obtain a sequence {— }^ 1 approaching 9 from below. In that 
case {n r 9} — > , and we get the estimate 

F + {k nr ) < n -{l + e){n r 9} < 7r(A + - ) 



2 v ' 1 r J 2^/1 



71,. 
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for any e > and r large enough. In the same way one can approximate 9~ l . By 
( |3.12| ) we find that there are infinitely many gaps if c > 5~ 1 / 2 ti 2 L~ 1 (1 + e) for any 
e > , i.e., the assertion (c). 

Let finally 9 be a Last admissible number. Without loss of generality we may 
suppose that 

limsup a 2n = oo 

holds for its continued-fraction representation; otherwise we use instead 

[0, a ,ax,...} ... a ^ 

[ai, a 2 , a 3 , . . .] ... ao = 







-i 



By ( |3.17|) there is a sequence {n r }^ 1 such that n^n,^} — > 0+ . This means that 
for all large enough r we have 

n r F + {k nr ) = n r tan ^— {n r #}^ < 7rn r {n r 6>} — > 0. 

Hence k nr F + {k nr ) — » too, so there are infinitely many values of accumulating 
at infinity for which kF + (k) < | . | 

Remark 3.3 The critical value of 7r 2 /v / 5 = 4.414... in the part (c) cannot be 
pushed down — cf. [^7], Sec.II.l]. The bound is saturated for the golden mean, 
9 = | (1 + y/5) , and moreover, there are no gaps if \c\L is below the critical value 
in this case. 



3.3 S' s lattice spectra 

Let us turn to lattices with a 5' s coupling. The relation (|3.6| ) can be rewritten as 

Dk 

— = —F(k; —vt, —v 2 ) , (3.18) 

so the bands are now determined by the inequalities 

Dk 

TF T (k) > ± — , ±D>0. (3.19) 

This concerns positive energies; on the negative halfline we have instead the condi- 
tion 

> ^ > (3.20) 

with -F±(k) given by (|3.13|) ; the change in sign is due to the fact that k is now 
in numerator on the left side of ( |3.18 ). Replacing F±(k) by similar expressions 
containing a single term, we obtain in this way the band spectrum of the one- 
dimensional array of 5' s interactions. 
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Remark 3.4 Notice that the latter coincides with that of the 8' array of the same 
parameters || Sec. III. 3], because the corresponding transfer matrices differ just by 
sign; in higher dimensions the relation between 8' and 8' s is not that simple. 

We shall use again the definitions ( p,14|) and employ the symbol a(£, 9, D) for 
the spectrum of the 8' s lattice Hamiltonian. The conditions (|3.19| ) and (|3.20| ) have 
the following easy consequences: 

Proposition 3.5 (a) The spectrum has a band structure. For any nonzero D the 
number of gaps is infinite, o~(£, 9, D) = U^Ljav, Pr] with a r < f3 r < a r +i . 

(b) a(£,9,0) = [0,oo) . 

(c) If D > , each a r equals Oj^ for A = i\ or £ 2 , and some n e Z . On the 
other hand, (3 r = ( 2I f) 2 with m E 2Z for D < and r > 2 . 

(d) ax = for D > , while «i < if D < . 

(e) For —£\ — £2 < D < we have (3\ < and a 2 = . 

(f) a(£, 9, D') n M + C a(£,9,D)f]R + holds for \D'\ > \D\ . 

(g) Each gap is contained in the intersection of a pair of 5' Kronig-Penney gaps - 
see Remark \3.^ — with the coupling constant D and spacings £\ and £2 ■ 

In distinction to the 8 lattice there are therefore always infinitely many gaps for 
D ^ . As in the Kronig-Penney case and its 8' analogue, the roles of bands and 
gaps are, roughly speaking, reversed. Comparing to the part (h) of Proposition |3.1| . 
however, the asymptotics of band widths is slightly more complicated. 

It is clear that if a band [oc r , j3 r ] with r ^> 1 is well separated from the rest of the 
spectrum, its width has the same leading term as in the Kronig-Penney situation, 

ft - Or = + 0{r~ l ) , (3.21) 



where £j is the length to which this band corresponds by Proposition |3.1| c. If 9 is 
rational and the points k n and k m coincide, we get in the same way 

(5 r - OL r = -| (£i 1 +£ 2 1 ) + Oir- 1 ) . (3.22) 

It may happen, however, that k n and k m are not identical but close to each other, 
so that they still produce a single band. It is obvious from ( |3.19| ) that this leads to 
an enhancement of the band width above the value given by (|3.22|) . The effect is 
most profound just before the band splits. Suppose, e.g., that D > and k m > k n 
with 8 ~. — k m — k n ~ 8/D£i . The band width in the momentum variable is then 
8+r] up to error terms, where 77 solves the equation 

rf - 89- l r) - 8 2 9- 1 = 0. 
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The enhancement due to band conspiracy is therefore 



S + V , m 29 + 1 + VT+49 

= 9{6) ■= ntn , ,x • (3-23) 



5(1 + 9 



The same can be done if the order of k n , k m is reversed; in combination with ( |3.21 ) 
and (|3.22|) this yields the following estimates on the band width, 

^ + Oir- 1 ) < (3 r -a r < 1 (£r 1 +^2 1 ) e(9) + O^ 1 ) , (3.24) 

where 

e{9) := max{g(9),g(9- 1 )}. 

It is easy to see that e{9) > 1 with \im g ^ e(9) = lim^oo e(9) = 1 . For 9 = 1 we 
have e{9) = j (3+^) = 1.309 . . . , but the strongest conspiracy occurs if the lattice 
spacing is close to two to one, and the wider of the conspiring bands is below the 
other one (above for D < ), because the right side of (|3.24) reaches its maximum 



at g(2) = | . The value 9 = 2 itself is, of course, integer so there is no enhancement. 
Summing this discussion, we have 

Proposition 3.6 The band widths of a 8' s lattice with D^O satisfy the asymptotic 
bounds ( \3.24 ), where e{9) < | . 



4 Conclusions 

We have said in the introduction that a choice of the coupling at graph nodes in a 
realistic model should follow from a suitable "zero-diameter" limit of a tube system 
of the same topology. To illustrate some problems which may arise, recall that a 
cross-type region in the plane with the Dirichlet boundary exhibits a bound state 



28 1 . Moreover, using the Dirichlet bracketing |25|, Sec. XIII. 15] in combination with 



the results of [|J, one can check easily that any branched (star-shaped) system of 



infinitely long tubes with the Dirichlet boundary has at least one bound state [17|; 
sometimes the number of bound states may be even large, e.g., for systems of a 
narrow-fork form as can be seen from 0]. This conclusion extends to situations 
where the connecting-region boundary is "rounded" (no squeezing allowed) and the 
tubes involved are only asymptotically straight |12|, [I7|, |18[ . 

The existence of a single bound state could be preserved in the zero-diameter 
limit provided the corresponding "coupling constant" is chosen with a proper sign 
as we have pointed out in Section 2.2. However, the couplings discussed here cannot 
accomodate multiple bound states. Using the more general boundary conditions 
( |1.2j ) does not help: it is easy to see that such a coupling has at most two bound 
states, €j = — with 

Kl '-=-A^B> " 2 ]= - A + (Ll)B ] 
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this happens if the denominators are negative. In general therefore a zero-diameter 
limit should be expected to work in an energy interval around the continuous spec- 
trum threshold which is kept fixed, or at most it remains small with respect to the 
intermode distances when the junction region is scaled. 

The analysis of both a single junction and a rectangular lattice shows the ex- 
ceptional role of the "free" operators, i.e., those having the 8 coupling with c = 
(which is the same as 8' with C = ), or 5' s coupling with D = . Their S-matrix 
elements, band profile on a lattice, and other properties are, of course, nontrivial due 
the branching; however, they are the simplest possible. For instance, the reflection 
and transmission amplitudes through an n wire junction are 

r = ±?—^, t = ±- (4.2) 

n n 

for the "free" 5 and 5' s coupling, respectively, independently of energy. Numerical 
calculations of transport properties through Y-j unctions [^|] suggest that at least 
for some systems of n coupled straight tubes, this might be the correct low-energy 
scattering limit. 

On the other hand, a junction of finite-width tubes can have various geometries; 
in fact, an experimentalist would hardly guarantee that three joined quantum wires 
have a perfect Y shape. Moreover, if the connection region supports a potential, 
albeit a weak one, the low-energy scattering properties would be substantially al- 
tered. Hence the "non-free" boundary conditions are also of physical interest; one 
can even conceive easily a tube system in which junction parameters are tuned by 
application of an external field. 

Up to now we had in mind mostly simple junctions. Replacing them by regions 
of a nontrivial topology we arrive at a situation to which the considerations of Sec- 
tion 2.3 might be regarded as a simplified picture. Since we have shown there that 
the 5' and S' s couplings are, at least within a fixed interval of "intermediate" ener- 
gies, modelled by complicated enough geometric scatterers, also the conditions fl2.5| ) 
and Q2.6|) are likely to have something in common with the real world. Moreover, we 
have seen that the coupling constant of the limiting ideal scatterer is nonzero and 
it is fully determined by the geometrical properties of its approximants. 

In this respect, a comment is due. Without giving any details, the authors of 



19] suggested that such "composed" junctions can be described by the boundary 



conditions (jl.lp with the "renormalized" parameter dependent on energy. This 
does not contradict our conclusions. For instance, in the simplest case n = 2 the 
corresponding reflection and transmission amplitudes, which differ just by the sign 
from ( |2.3|) , may be written formally as the corresponding 8 scattering quantities 
provided we choose c(k) := —Dk 2 . Hence if one wants to describe a "composed" 
junction by means of a "dressed" coupling constant — which anyhow makes sense 
only when a prescription to compute the latter is given — it may happen that it 
differs substantially from the "bare" coupling. 

With this we leave this subject and turn to latice Hamiltonians of Section 3. 
In addition to their possible use as models of quantum wire superlattices, they 
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represent an interesting mathematical object, and the observed dependence of the 
spectra on number-theoretical properties of the parameter 9 raises many questions. 
One would like to know, for example, how the band and gap patterns do actually 
look, what are their fractal properties with respect to 9 , or what is the measure of 
the spectrum relative to a suitable measure on M. We intend to return to these 
problems in a later publication. 

The results of Sections 3.2 and 3.3 show that despite there are "less" gaps in two 
dimensions, and despite the behaviour of bands and gaps is in general irregular — for 
an irrational 9 they exhibit asymptotically a "squared quasiperiodic" distribution 
- it coincides roughly with that of the Kronig-Penney analogues to our lattices, 
namely that for the 5 coupling the bands dominate at high energies, while the 
converse is true for the S' s . In analogy with the one-dimensional case [§, [7], |23[ , 
one can therefore make a conjecture concerning the situation when a 5' s lattice with 
D 7^ is placed into an electric field. The heuristic tilted-band picture suggests 
the existence of localization; the spectrum will remain continuous, of course, but an 
unrestricted propagation may be possible only in the direction perpendicular to the 
electric field. In the 5 lattice case, where we have for the one-dimensional situation 
a guess but no rigorous result, the problem is even more exciting; the results of the 
present paper show that at least for some values of the lattice parameters there is 
no localization. 
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Figure captions 

Figure 1. Scattering on T n (N,£) with n = 3 and N = 2 . 
Figure 2. A rectangular latice. 
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This figure "figl-l.png" is available in "png" format from: 
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